wo 2005/043138 



- 1 - 



PCT/IB2004/003559 



A TIME-DOMAIN METHOD AND APPARATUS FOR 
DETERMINING THE DEPTH AND CONCENTRATION OF A 
FLUOROPHORE IN A TURBID MEDIUM 

5 TECHNICAL FIELD 

This application relates to fluorescence 
measurements in turbid media. 

BACKGROUND OF THE INVENTION 

The use of fluorophores for identifying or 
localizing structures within objects is of considerable 
importance in diverse fields and in biological/medical 
applications in particular. Optical fluorescence imaging 
of turbid media such as biological tissue, is primarily 
achieved with Continuous Wave (CW) methods. Typically a 
light source is employed to illuminate the object of 
interest^ e.g. a mouse in vivo^ and the emitted CW 
fluorescence intensity signal from the flurophore is 
measured directly with a camera. Systems in which the 
fluorescence from small animals comprising bioluminescent 
molecules can be measured have been described. See for 
example, U.S. Patents 6,775,567 and 5,650,135. 

For applications such as diagnosis and 
pharmacological studies, it is often desirable to 
25 determine the concentration and depth of the fluorophore. 
However, to assume that the direct CW fluorescence 
intensity signal is proportional to the flurophore 
concentration can be misleading since the depth of the 
flurophore will also impact the CW fluorescence intensity 
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signal. In fact, given a single CW source and CW detector 
measurement - it is impossible to decouple f lurophore 
concentration and depth. 

5 To address this problem of decoupling depth and 

concentration in CW measurements, some researchers are 
exploring CW tomography, in order to reconstruct a three- 
dimensional distribution of the flurophore concentration 
in the object, thereby attempting to account for the 

10 depth of the flurophore. However, these tomographic 
approaches require multiple source-detector pair 
measurements from many angles, combined with complex, 
computer intensive inversion algorithms. Furthermore, CW 
tomography requires an assumption about the scattering 

15 coefficient of the object, since CW can neither decouple 
the intrinsic absorption and scattering of the object. 

The most complete description of photon migration in 
turbid media is provided by Time Domain (TD) optical 

20 methods which have previously been used to decouple the 
attenuation coefficient, given from CW ihtensity 
measurements, into the underlying absorption and 
scattering coefficients-. However TD methods have not been 
applied to obtain depth and concentration of 

25 f luorophores . 

SUMMZ^Y OF THE INVENTION 

The invention relates to the measurement of 
fluorophores localization and concentration within a 
turbid medium. More specifically the invention relates to 
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the measurement of fluorophores depth and concentration 
within a turbid medium. 

There is advantageously provided a method wherein TD 
optical methods are used to decouple fluorophore depth 
5 and concentration from a single source and detector 
measurement employing a direct analysis technicjue^ 
thereby avoiding the aforementioned burden of CW 
tomography and overcoming the limitations of conventional 
CW direct fluorescence intensity imaging. 

10 In one embodiment there is provided method for 

determining the depth of a volume comprising a 
fluorophore in a turbid medium using time domain (TD) 
optical fluorescence by obtaining Temporal Point Spread 
Function (TPSF) data, by injecting light at an injection 

15 point on an object at an excitation wavelength of the 
fluorophore and detecting light at a detection point of 
the object at an emission wavelength of the fluorophore, 
determining a time (tmax) substantially corresponding to 
the maximum of the TPSF and correlating tmax with the 

20 depth. In a preferred embodiment, the injection and 
detection points are a back- reflection geometry and are 
substantially equidistant from the volume comprising the 
fluorophore. In this case, the depth may be expressed* as 
a function of the speed of light in the medium, the 

25 scatter coefficient in the medium, the lifetime of the 
fluorophore in the medium and tmax- Thus, by measuring or 
providing an estimate of the above parameters, the depth 
of fluorophore can be estimated. 
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In another embodiment there is also provided a 
method for estimating the concentration of a fluorophore 
in a volume of a turbid medium using optical fluorescence 
by obtaining the depth of the volume comprising the 
5 fluorophore, providing optical properties for the medium, 
obtaining a intensity surface reflection measurement of 
the fluorophore, and normalizing the CW intensity 
measurements using optical properties of the medium and 
the fluorophore depth to obtain a relative fluorophore 
10 concentration. An absolute concentration may also be 
obtained by calibrating the fluorescence emission using 
known concentrations of the fluorophores embedded in the 
medium. 

In yet a further embodiment there is provided a 
15 method for generating a tomographic image of a 
fluorophore distribution in a turbid medium by obtaining 
a topographic image of the fluorophore distribution, 
determining the depth of a plurality of volumes of 
interest comprising the fluorophore using depth 
20 determination methods of the present invention, and 
combining the depth information and the topographic image 
to generate a tomographic image of the fluorophore 
•distribution . 

In another aspect of the invention there is provided 
25 an apparatus for determining the concentration of a 
fluorophore in a turbid medium which comprises a light 
source, optically coupled to a source channel and the 
object, to inject light in said object at a desired point 
and excitation wavelength; a first detector channel, 
30 optically coupled to a photon detector and the object, in 
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a back-reflection geometry relative to the source 
channel, to acquire at least one TPSF from a desired 
point on the object to determine the depth of the 
fluorophore, a second detector channel in a trans- 
5 illumination geometry relative to the source channel, to 
measure an emission intensity of the fluorophore, a means 
for spatially positioning the object relative to the 
channels, a depth calculator and a concentration 
calculator. 



10 BRIEF DESCRIPTION OF THE DRAWINGS 

Further features and advantages of the present 
invention will become apparent from the following 
detailed description, taken in combination with the 
appended drawings, in which: 

15 Fig. 1 is a 3-dimensional graphic relating tmax/ 

depth of fluorophore and the concentration of the 
fluorophore; 

Fig* 2 is a 3-dimensional graphic relating CW signal- 
intensity,, depth of fluorophore and the concentration of 
'20 the fluorophore; 

Fig. 3 is a schematic representation of a preferred 
fluorescence signal acquisition configuration for depth 
determination; 

Fig. 4A is a schematic representation of the back- 
25 reflection geometry; 

Fig. 4B is a schematic representation of the trans- 
illumination, geometry; 
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Fig. 5 is a schematic representation of an example 
of the system combining the back-reflection and the 
trans-illumination geometries; 

Fig. 6 is a schematic representation of an example 
5 of the apparatus combining the back-reflection and the 
trans-illumination geometries; and 

Fig. 7 is a schematic representation of an example 
of the apparatus combining the back-reflection and trans- 
illumination geometries wherein the trans-illumination 
10 detector channel is a CCD camera. 

It will be noted that throughout the appended 
drawings, like features are identified by like reference 
numerals. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
15 The fluorescence emission from a fluorophore 

embedded in a turbid medium can be measured over time 
using time domain (TD) optical fluorescence to generate 
Temporal Point Spread Function - (TPSF) data. Acquisition 
of such data will be described in more details below^ but 
20 briefly it consists in injecting a pulse of light at an 
excitation wavelength in the medium comprising the 
fluorophore and detecting the emission fluorescence as a 
function of time. TD fluorescence experiments were 
conducted with a single source-detector measurement for a 
25 small inclusion. of fixed flurophore concentration 
submerged into a turbid, liquid medium. It was observed 
that the temporal position, of the maximum of the TPSF 
generated by the fluorescent photons (tmax) increased with 
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the depth of the submersion of the inclusion in the 
turbid medium. 

TD fluorescence experiments were also conducted with 
single source-detector measurements for small inclusions 
5 of various fluorophore concentrations' submerged at fixed 
depths in a turbid, liquid medium. It was observed that 
tmax generated by the fluorescent photons remained constant 
for inclusions of the same submersion depth, regardless 
of their fluorophore concentration. These results are 

10 displayed in figure 1 which is a tri-dimensional plot of 
tmax as a function of depth and concentration of the 
fluorophore. Hence it was demonstrated that the 
measurement of tmax from a fluorophore inclusion allows 
the determination of the depth of a fluorophore 

15 independent of its concentration. 

For the above series of experiments, CW data was 
also generated by temporally integrating the complete 
TPSF. As . expected, the CW fluorescence intensity 
decreased as the depth of the fixed fluorophore 

20 concentration increased, but it also decreased as the 
fluorophore concentration of inclusions at a fixed depth 
decreased, as shown in Figure 2. Hence, it was 
demonstrated that direct CW fluorescence intensity 
measurements alone cannot decouple fluorophore inclusion 

25 concentration from depth. 

In the present description it will be appreciated 
that the intensity of the fluorescence emission of the 
fluorophore can be obtained using direct CW measurements 
or by integrating a Temporal Point Spread Function (TPSF) 
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obtained by the time domain (TD) method. It will be noted 
that time domain data can be obtained by connecting 
frequency domain data using Fourier Transform. 

It will also be appreciated that the term 
5 fluorophore can refer to either an extrinsic fluorophore 
which is understood to be a fluorophore that is added to 
a given medium or an intrinsic fluorophore which is 
understood to be a fluorophore that is normally comprised 
in a given medium. For example biological tissue may 
10 comprise molecules that naturally fluoresce and are 
therefore intrinsic f luorophores . 

In the present description it is. also meant by 
turbid medium, any medium that scatters light. Such media 
may comprise but are not limited to biological tissues 
15 such as breast,- brain, muscles and the like. 
Furthermore, while the embodiments described herein 
examplify the application of the- method and apparatus to 
. small animals, the method and apparatus can also be 
applied to larger animals including humans. 

20 Thus in view of the results shown in figure 1 and 2, 

the depth, within an object comprising turbid medium, of 
a volume comprising a fluorophore can be determined by 
establishing a calibration curve for which tmax 
measurements are obtained for known depth of the volume 

25 and using this calibration curve to estimate the depth of 
a volume comprising a fluorophore which depth is 
originally unknown. 

The depth of more than one fluorophore can be 
assessed by injecting light of different wavelengths 
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corresponding to the excitation wavelengths of the 
different fluorophores and detecting light at an emission 
wavelength for each of the fluorophore. 

Under certain assumptions such as assimiing that the 
5 optical properties of the medium are the same at the 
excitation and emission wavelength^ the fluorescence 
intensity as a function of time can be expressed by the 
Born approximation: 



r _±\ 

e ^ 



\ J 



*(ZR) (equation 1) 



10 Where: 

rgp is the distance from source s (point on the 
object at which light is injected) to fluorophore depth 
position p; 

rpd is the distance from fluorophore depth position p 
15 to detector d; 

|ia is the optical absorption coefficient; 



D is the optical diffusion coefficient, £) = where 

; |is' is the reduced optical scatter coefficient; 
V is the speed of light in the medium; 



20 



Q is the quantum efficiency; 



C is. the concentration of the fluorophore; 
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T is the lifetime of the fluorophore; and 

IR is the impulse response of the instrument used to 
measure fluorescence 

By setting the first derivative, of equation 1 as a 
5 function of time equal to zero, the time position of the 
maximum of the TPSF (tmax) can be found. Under certain 
approximations (absorption is small at time shorter than 
tmaxf the scatter coefficient is known or can be 
• approximated) and by assuming that rgp is approximately 
10 equal to rpd, it is found that the following equation can . 
be derived from equation 1: 

t max = ^^^-^ — (equation 2) 

where d is the depth of the volume comprising the 
fluorophore . 

15 Figure 3 schematically depicts how the TPSF may be 

acquired. One of more pulse of light at an excitation 
wavelength is injected at injection point 10 in an object 
(in this case examplified by a mouse) comprising one or 
more fluorophore. The photons diffuse within the tissue 

20 and a certain fraction of the photons eventually reach 
the fluorophore molecules contained in a volume of 
■ interest (VOX) 12 which are thereby excited. The 
molecules then emit a fluorescence signal. The emission 
photons will diffuse and a fraction of them will reach 

25 detection point ' 14 to produce TPSF data 16. In a 
preferred embodiment the acquisition of the fluorescence 
signal is performed in the back-reflection mode as shown 
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in figure 3, with the injection and detection points 
substantially equidistant from the volume comprising the 
fluorophore. That is rgp 17 and rpd 18 are approximately 
equal • 

5 As can be seen from equation 2, the depth of the 

fluorophore is proportional to tmax but also to the speed 
of light, the scatter coefficient and the lifetime of the 
fluorophore. In order to deteimiine the depth of a volume 
comprising a fluorophore in a turbid medium, the speed of 
10 light, scatter coefficient and lifetime can be provided 
using known values for the medium. Alternatively, the 
scattering coefficient and lifetime can be provided using 
the standard time domain approach by direct measurements. 

The scattering coefficient can be determined at 
15 either the excitation and/br emission wavelengths of the 
fluorophore, since in a preferred embodiment the method 
assumes that these values are the same for the excitation 
and emission wavelengths. However, it is also possible to 
derive tmax by using the scattering coefficient as 
20 determined at both the excitation and emission 
wavelengths of the fluorophore. In this respect, it will 
be appreciated that a more general form of equation 1 in 
which the absorption and scatter coefficients, determined 
at both the emission and excitation wavelengths, are used 
25 is also encompassed in the present invention. 

The calculation from which equation 2 is derived is 
based in part on the assumption that the distance between 
the point of light injection 10 at the surface of the 
object and volume 12 and between volume 12 . and the point 
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of fluorescence emission detection 14 is substantially 
identical. Determination of the injection and detection 
points relative to the volume of interest in the 
reflection geometry can be achieved, for example by 
5 obtaining a topographic image of the fluorophore in a 
region of interest. While the topographic image can be 
acquired using the system of the present invention,- it 
can also be acquired by systems and methods that are well 
known in the art (see for example WO 2004/044562 and US 
10 patent 6,217,847) . In this respect the 2D (topographic) 
image may be obtained by optical modalities such as TD, 
CW and frequency domain. (FD) . 

The injection and detection points may also be 
determined based on prior knowledge of the expected 
15 location of the fluorophore in a plane substantially 
perpendicular to the depth. For example, the location of 
the fluorophore may be known and this information may be 
used to select the injection and detection points. 

Once the injection and detection points have been 
20 determined, TPSF data is then acquired by injecting 
pulses of light at an excitation wavelength of the 
fluorophore at the injection point and by detecting as a 
function of time the fluorescence emission at an emission 
wavelength. The time at which the TPSF reaches a maximum 
25 (tniax)is then recorded and used in the above relationships 
(either equation 1 or equation 2) to determine the depth 
of the fluorophore. 

It ' will be appreciated that a single 
injection/detection points pair is sufficient to 
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determine the depths therefore providing a rapid method 
for determination of depth- However, it is also possible 
to use multiple injection/detection pairs to improve the 
accuracy of depth determination. 



signal detected at the surface can be related to 
fluorophore concentration by the optical properties of 
the medium (absorption and scattering coefficients) and 
the depth of the fluorophore. 



Using time-domain information as described above, 
the depth d can be determined. Isolating C in equation 3 
and knowing signal intensity and depth of the 

15 fluorophore, one can thus recover the concentration of 
fluorophore (i.e. the amount of fluorescent molecules per 
unit volume) within an accuracy that depends 
exponentially on the recovered depth accuracy- Thus, in 
another aspect of the invention, estimates of the 

20 relative concentration of the fluorophore. Cone . Relative/ 
can be obtained by determining .its depth, d, and 
normalizing the surface intensity measurement, J, as 
follows (Equation 4) : 



5 



For a given depth, the intensity. J of the emission 



10 




(equation 3) 



Cone, 



Violative 



(equation 4) 



25 



under certain assumptions. Equation 4 can be derived 
from equation 1 . 
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Determination of the concentration may also be 
derived from equation 1. Referring to figure 4a, it can 
be seen that, when measurements are performed in a back- 
reflection geometry, the distances rsp 17 and rpd 18 depend 
5 on the depth d of the fluorophore. Since there is a 
certain variation associated with the measurement of the 
depth determination, the determination of the 
concentration will carry the variation associated with 
depth determination. Therefore, in a preferred 
10 embodiment, the intensity J may be acquired in the trans- 
illumination mode as shown in figure 4B, whereby the 
distance between the point of injection 10 and the point 
of detection 15 is constant for a given object and is 
equal to rgp + rpd. Thus, in the trans-illumination mode 
15 the term ^^rsp + rpd" does not depend on d. It is therefore 
possible to obtain a better estimation of the 
concentration . 

In a further embodiment of the invention there is 
therefore provided a method for determining the 
concentration of a fluorophore comprising obtaining the 
depth of the fluorophore using a back-reflection geometry 
as described above and a fluorescence intensity 
measurement in a trans-illumination geometry and 
calculating the concentration based on the measured 
depth . 

This method provides a means for rapidly comparing 
the relative concentration of a fluorophore at two or 
more distinct locations. The determination of the 
30 relative concentration of a , f luorescently labeled 
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antibody in different organs of an animal provides a non- 
limiting example. 

Relative concentration can be further processed, to 
estimates of absolute concentration if the surface 
5 intensity measurement is absolutely calibrated. That is 
to say a calibration curve can be established for J using 
known concentrations of the fluorophore. The scattering 
and absorption coefficients can be provided using known 
values for the medium or alternatively may be provided 

10 using the standard time domain approach by direct 
measurements. Biological tissues may not conveniently 
allow for easy calibration since incorporating known 
concentrations of fluorophores may be difficult. In these 
cases the calibration may be effected in a turbid medium 

15 having similar optical properties as the biological 
tissue to be measured. 

By further analyzing the fluorescence intensity data 
from the above experiments (figure 1 and 2)^ it was 

20 possible to calibrate the optical attenuation of the 
turbid medium. In other words, for the inclusion of fixed 
( fluorophore concentration, the optical attenuation can be 
obtained by fitting the decrease of fluorescence 
intensity measured by the surface detector as the 

25 inclusion depth increases. The same optical attenuation 
can then be used to describe the decrease of fluorescence 
intensity of fluorophore with other concentrations as the 
inclusion depth increases. 

Given the above series of experiment, the inclusion 
30 may be submerged at unknown depth and concentration in 
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the turbid medium. As stated above, direct fluorescence 
intensity measurements alone cannot decouple the 
inclusion depth from its fluorophore concentration. 
However TD experiments can provide tmax which yields the 
5 inclusion depth. Once the inclusion depth is known, the 
calibrated optical attenuation of the turbid medium can 
be used to normalize the fluorescence intensity 
measurement to yield the relative concentration of the 
fluorophore inclusion. Further, an absolute calibration 
10 of the intensity measurements permits estimates of the 
absolute concentration. 

In yet a further embodiment of the invention there 
is provided a method for generating a tomographic image 
of a fluorophore distribution. Topographic images of a 

15 region of interest can be obtained using well known 
method such as CW, TD frequency domain (FD) optical 
imaging. Such topographic images can be combined with the 
depth and/or concentration information obtained by the 
methods described above to generate a 3-dimension 

20 (tomographic) image of a volume of interest comprising 
one or more f luorophores . 

It will also be appreciated that the surface 
intensity may be provided by temporarily integrating the 
25 time domain fluorescent measurement. 

With reference to figure 5, an example of the. system 
of the present invention for measuring fluorophore 
concentration is schematically illustrated. Light source 
50 is optically coupled to object 51, which comprises the 
30 fluorophore, through source channel 52. Light emitted by 
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the fluorophore is optically coupled to photons detector 
53 through back-reflection channel 54. When operated in 
the TD. mode, the measured TPSF can be acquired and used 
by depth calculator 55 to determine the depth of the 
5 fluorophore. The intensity of the fluorescence signal can 
be measured either simultaneously with the TPSF or 
sequentially by detecting photons with photons detector 
56 optically coupled to the object through trans^ 
illumination channel 57. The measured intensity is then 

10 combined with the calculated depth to calculate the 
concentration at concentration calculator 58. Calculation 
of the depth can be performed with equation 2 which 
requires estimation of scatter coefficient, speed of 
light in medium, and lifetime of fluorophore which are 

15 provided by estimator 59. Acquisition controller 60 
controls source/detector geometry and source 
characteristics such as intensity, pulse length, 
wavelengths and the like. 

A schematic diagram of an embodiment of the 
20 apparatus is shown in figure 6. Source 61 provides light 
at the excitation wavelength of a fluorophore embedded in 
object 62. The light is directed towards a predetermined 
point of light injection on object 62 using source 
channel 64. The source channel is an optical means for 
25 directing the excitation light to the desired point on 
the object and may consist of a fiber optic, reflective 
mirrors and the like. A first detector channel 65 is 
positioned to detect emission light in a back-ref lection 
geometry and a second detector channel 66 is position in 
30 a trans-illumination geometry. The detector channels are 
optical means for collecting the emission light from 
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desired points on object 62 and are optically coupled to 
photon detection 69. As for the source channel, the 
detector channels may consist of a fiber optic, 
reflective mirrors and the like . The source and detector 
5 channels may operate in a contact or free space optic 
configuration. By contact configuration it is meant that 
one or more of the components of the source and/or 
detector channel is in contact with object 62. By free 
space optic it is meant that light is propagated through 

10 air and directed to or collected from the desired points 
with appropriate optic components. An example of a free 
space optic system is described in WO 2004/044562. The 
detector channels can be coupled to spectral filters 67 
to selectively detect an emission wavelength of . the 

1 5 f luorophore . 

The source and detector channels can be physically 
mounted on a common gantry 68 so as to maintain them in a 
fixed relative position. In such an arrangement, the 
position of the ppint of injection of excitation light 

20 and that of the point from which the emission light is 
collected can be selected by moving (scanning) the gantry 
relative to the object. Alternatively relative 
positioning of the object and the source/detector 
channels may be accomplished by moving the object 

25 relative to the gantry. 

The position of the source and detector channels may 
also be controlled independently from one another. It 
will be appreciated that the position of the back- 
reflection and trans-illumination detector channels can 
30 also be independently controlled. Furthermore the system 
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may also allow a combination of arrangements . For 
example, the trans-illumination channel may be in a fixed 
position relative to the source channel whereas the 
position of the back reflection channel is controlled 
5 independently. The favored arrangement may depend on the 
type of object being probed, the nature and/or 
distribution of the f luorophore and the like . In a 
preferred embodiment, the source channel and detector 
channel are co-axially positioned. 

10 While the back-reflection detection is preferably 

done in the time-domain to provide depth information, the 
intensity of the fluorescence in the trans-illumination 
geometry may be detected by either CW, TD or FD domain 
modality. 

Detectors for the emitted light in the time domain- 
may comprise time correlated single photon counting 
devices (TCSPC) , time gated intensified charged couple 
device (ICCD) and the like. For continuous ' wave 
measurements photon multipliers, CCD cameras and the like 
may be used. The detector channel in the trans- 
illumination geometry may comprise a CCD camera 71 as 
show in figure 7 . 

The object can be placed on a transparent platform 
70 or can be suspended in the desired orientation by 
25 providing attachment means that can be removeably secured 
to the object and an appropriate structure within the 
apparatus. Preferably the position of the platform or the 
attachment means can be adjusted along all three spatial 
coordinates . 
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In the trans-illximination geometry the thickness of 
the object is preferably determined to provide a value 
for the optical path (rgp + rpd) . If the channels are in a 
contact configuration the thickness may be provided by 
5 the distance between the source channel at the point of 
light injection and the trans-illumination detector 
channel at the point of light collection. In the case of 
a free space optic configuration . the thickness may 
advantageously be provided by a profilometer which can 
10 accurately determine the coordinates of the contour of 
the object. 

A plurality of fluorophores can be monitored 
simultaneously by providing appropriate excitation and 
emission wavelengths. To this end source 61 may consist 

15 of a plurality of sources operating at different 
wavelengths. Alternatively the source may be a broadband 
source optically coupled to a spectral filter to select 
appropriate wavelengths. The emission wavelengths can be 
de-multiplexed into individual wavelengths by spectral 

20 filter 69. Selection of wavelength may also be effected 
using other appropriate optical components such as 
prisons . 

As described above the depth of the fluorophore may 
be determined by measuring tmax from a TPSF obtained from 

25 time-domain measurements. Furthermore the concentration 
of the fluorophore is advantageously measured using a 
trans-illumination geometry. In a preferred embodiment, 
the apparatus of the present invention provides for the 
simultaneous measurement of the depth, using a back- 

30 reflection geometry, and the concentration using a trans- 
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illumination geometry. Advantageously, this simultaneous 
determination can be achieved with a single source of 
light, and the positioning of the source channel can be 
coordinated with the positioning of the back-reflection 
5 and the trans-illumination detectors. 

It will be appreciated that multiple source and 
detector channels can be provided for simultaneous or 
sequential multiple measurements. 

The embodiment (s) of the invention described above 
10 is (are) intended to be exemplary only. The scope of the 
invention is . therefore intended to be limited solely by 
the scope of the appended claims. 



